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Abstract 

We perform Monte-Carlo calculations of the radiative transfer of Lya photons emitted by a source embedded in 
a neutral collapsing gas cloud. This represents a young galaxy or quasar during the early stages of the epoch of 
reionisation (EoR). After computing the Lya spectrum as function of radius and time, we find that the Lya color 
temperature may be negative in large volumes surrounding the central source. Motivated by this result, we explore 
the prospects for a population inversion in the hyperfine levels of atomic hydrogen via the Wouthuysen-Field (WF) 
effect. The reason for this exploration is clear: if 21-cm masers exist during the EoR, they could greatly boost the 
expected 21-cm flux from this epoch. We find that population inversion is unlikely to occur for several reasons: (1) 
the required Lya luminosities of the central source exceed ~ 10^^ erg s~^. The radiation pressure exerted by such 
a large Lya flux likely halts the collapse of the cloud; (2) When quantum corrections to the WF-coupling strength 
are applied, the required Lya luminosities are (even) larger by orders of magnitude; (3) A relatively low flux of Lya 
photons that is produced via other channels (x-ray heating, coUisional excitation of hydrogen) prevents the Lya color 
temperature from becoming negative. 
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1. Introduction 

Detecting the redshifted HI 21-cm line from the 
epoch of reionisation (EoR) is among the greatest 
observational challenges in present-day cosmology. 
Experiments such as the Low Frequency Array (LO- 
FArII], the Primeva l Structure Telescope (PAST, 
IPeterson et all 20061 ) and the Murchinson Widefield 
Array are currently being developed with 
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the goal of detecting this EoR 21-cm signal within 
the next few years. 

The strength and sign of the 21-cm signal depend 
on the spin temperature of the hydrogen gas, T^, 
which is defined through the ratio of atoms in the 
ground and excited states of the hyperfine transition 
in the electronic ground state of hydrogen, n2/ni = 
52/31 exp[— r,/Ts], where the ground and excited 
states are denoted by the labels '1' and '2', respec- 
tively and = 0.068 K. If we denote the tempera- 
ture of the Cosmic Microwave Background (CMB) 
by TcMB , then neutral hydrogen gas can be seen in 
emission if Ts > Tcmb ■ When < Tcmb , the neu- 
tral hydrogen gas can be seen in absorption against 
the CMB. 

The distribution of atoms in the hyperfine levels 
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is affected by collisions between atoms, absorption 
and re-emission of the Cosmic Microwave Back- 
ground (CMB) photons, and absorption a nd re - 
emission of Lya photons ( Furlanetto et all l2006h . 
These Lya photons were either emitted by nearby 
sources, produced as local recombination emis- 
sion, or in cascades following photo-excitation 
by photons with energi e s E g [10.2 — 13.6] eV 



(jPritchard fc Furlanettol . l200d ). When the Lya 



more detail whether < barely needs justifica- 
tion: if high redshift 21-cm masers exist, they could 
boost the 21-c m EoR signal by orde r s of m agnitude. 
Interestingly, ICbizhov fc Shanirol ^oM) recently 
showed that scattering of Lya and Ly/? photons 
by deuterium can cause its spin temperature to be 
negative. Technically, this produces a 92-cm maser, 
although the associated gain is too small to provide 
a dramatic effect. 

Bef ore proceeding, we point out that jPeguchi fc WatsonI 
(jl985f ) found that Lya scattering in expanding gas 



scattering rate dominates the coUisional and CMB 
scattering rate, the hydrogen spin temperature, 

Tg, approach e s the Lya color temperature, Ta clouds results in ^ Tgas. One expects the same to 

apply to a collap sing gas cloud. However, t he anal- 
ysis presented in Deguchi fc Watson ( 1985[ ) focused 
on Lya that was generated throughout the cloud. 
On the other hand, the young pre-reionisation 
universe may contain bright Lya sources, such as 
young galaxies and/or (mini)-quasars, that are sur- 
rounded by cold, neutral, collapsing gas. In such 
scenarios the Lya emission is concentrated strongly 
towards the center of the collapsing cloud. This 
introduce s an important difference with the calcu- 
lations of iDeguchi fc Watson Photons that 
are emitted in the center of the cloud must propa- 
gate through a larger column of collapsing material. 
This causes more energy to be transferred from the 
infalling gas to the photons, which re sults in a larger 



(|Wouthuvsenl . Il952l ). The color temperature has 
been shown to approach the gas temperature, 
Ta Tgas, because of the frequency redistribution 
that Lya photons experience while resonantl y scat- 
tering through an optically thick medium ([Field, 
Because of this 'Wouthuysen-Field' effect, 
it is commonly assumed that the presence of Lya 
photons strengthens the coupling of the spin to the 
gas temperature (for the most detailed calculation 
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see 

detectable in emission. 

In this paper, we examine whether there are 
physical conditions under which the WF-effect pro- 
duces population inversion (i.e. Tg < 0) in the 
hyperfine levels. This investigation is motivated 
by recent Monte-Carlo calculations of Lya radia- 
tive transfer, which show that Lya photons emerge 
from neutral colla psing gas clouds with a system- 
atic blueshift ( e.g.lZhengfc Mir alda-Esc udi 12002 : 
Diikstra et all . l2006l ). This blueshift was found to 
increase with Ly a optical depth and w ith the infall 
speed of the gas (I Piikstra et al. , 20061 ) . In this pa- 
per we show that this may lead to a negative Lya 
color temperature, and potentially to a negative 
spin temperature. The possibility of population in- 
version via the WF-effect and the accompanying 
stimulated 21-cm e mission has been addressed in 
previ o us work (e.g. Var shalovich. '1967"; 'Shklovskii', 



overa ll blue-shift of the Lya line (jPiikstra et al 



1967 : Ivan Bueren fc Oo rt. 1968; Kaplan et aL, 



19791), but (obviously) not in the light of the re- 
sults of these recent radiative transfer calculations, 
and not in the context of the EoR Q- Exploring in 



20061 ). and we will show in this paper that this large 
blueshift is the main reason that the WF-effect may 
cause a population inversion. 

The outline of this paper is as follows. In §[5] we 
discuss the expected maser amplification associated 
with a given negative hydrogen spin temperature. In 
§ |3] we derive an expression for the Lya color tem- 
perature, which can be applied to arbitrary Lya fre- 
quency distributions, and its relation to the hydro- 
gen spin temperature. In §[4] we discuss an example 
of a Lya source surrounded by cold collapsing gas, in 
which the Lya color temperature is negative. This is 
followed in §[S]-§[n]by a discussion of the factors that 
determine whether 21-cm masers can actually exist 
inside such clouds, before presenting our conclusions 
in § [T] The parameters for the background cosmol- 
ogy used throughout this pap er are flm = 0.3, f^A 



0.7, nb = 0.044, and h = 0.7 (ISoergel et all . 120031) 



* ISpaana &: NormarJ l|l997l ) have investigated whether hy- 
drogen masers can exist during EoR. However, they focused 
on hydrogen na recombination lines, which correspond to 
the transitions of the form n — > n — 1, w here n is the princi- 
pal qu antum number of atomic hydrogen. ISpaans &6 NormanI 
lll997h found that maser lines may be produced for n ~ 120 



along the edges of early HII regions. 



2. The Gain Factor of a 21-cm Maser 

Stimulated emission of 21-cm emission occurs 
when Ts < 0. Here, we estimate the expected gain 
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in the total 21-cni emissivity. The change of the in- 
tensity in the hne cent er of the 21-cm hne, J, with 



distance s, is given by ijRvbicki fc Lightmanl [l979h 



/(s) = /(0) e-"^ + ^(l 
a 



(1) 



where a is the opacity of the gas at As =21cm and 
j is its volume emissivity. Furthermore, /(O) is the 
intensity of the 21-cm photons before they enter the 
gas, /(O) = 2fcsTcMB/A|. The 21-cm volume emis- 
sivity, J, and opacity, a, are given by 

J = -^7.2^21 (2) 
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where 1^5 = 1.4 Ghz, is the frequency of the hyper- 
fine transition; Avs — vsvth/c, in which wth is the 
thermal velocity dispersion of the atoms; A21 and 
B21 are the Einstein coefficients associated with the 
hyperfine transition. Eq. ([1]) can be rewritten as 

I[s) = /(0)e^/^- - (e^/«- - 1) , (3) 

where Rm = -l/a. We used B21 = A2ic'^ /2hvl, 
ni ^ ?T._ff/4 and the definition of the spin tempera- 
ture (§[1]) to rewrite j/a = 2kBTs/X^. Furthermore, 
we used that <ti \Ts\, which is valid throughout 
this paper. Since for a negative spin temperature 
a < 0, /(s) increases exponentially with s. Rm is the 
length scale over which I{s) increases by one factor 
of e and can be written as 

z\-3/l + ,5\-i 



(4) 



Thus a reasonably small region in which T, < 0, 
also known as the "masing region" , could produce 
a significant boost in /(O). According to Eq. (g]), 
the total maser amplification increases toward lower 
Tgas and \Ts\ values and higher densities. 

3. The Lya Color and Hydrogen Spin 
Temperature 

3.1. The Lya Color Temperature 



as 



Madau et al.l (jl997f ) defined the color temperature 

1 dJaii^)- 





raiogACl 


-1 r 









Ja.(y) dv 



where Mu ~ c^Ja {v)/'2,hv^ is the photon occupation 
number. The more general expression for the Lya 
color temperature is f Appendix I A. ip 

kTg _ ^ J Ja{v)4>{v)dv _ J Ja{v)(l){v)dv 



(6) 



(also see Meiksin, 2006, Eq. 16). Here, (t){v) is 
the normalised line profile, / <l)[u)dv — 1 (e.g. 
iRvbicki fc Lightniaii . 19791 ). Frequency is expressed 
in terms of a dimensionless frequency variable, 
X = {v — vq)/ AvD, where Aud = vovth/c, in which 
j/Q — 2.46 X 10^^ Hz is the Lya frequency. A nega- 
tive/positive value of X corresponds to a Lya pho- 
ton that is redshifted/blueshifted with respect to 
the line center, respectively. In the 'core' of the line, 
(j){x) ~ , while in the 'wing' 0(a;) ~ a/[^x^]. 
Here, a is the Voigt parameter, which is given by 
a = Aa/AnAiyD = 4.7 x lO""^ (13 km s" V-fth), 
where Aa is the Lya Einstein-A coefficient. The 
transition between core and wing is defined to occur 
at the frequency at which e~^* = a/[-</TrXf]. Since 
the function (t>{x) is strongly peaked around the 
line center, photons in this frequency range usually 
dominate the contribution to the integral. 

At first glance the color temperature becomes neg- 
ative when there are more photons with ly > vq 
than with v < vqi because d(j)/du < for v > vq. 
This may occur in gas clouds moving towards a Lya 
source, as these w ould 'see' a blueshifted Lya line 
(jShklovskiil . Il967r ) . As was mentioned in § [T] how- 
ever, repeated scattering of photons in the line core 
rearranges the m in frequency such tha t Ta — *■ T^ 



IField . 1959 ; van Bueren fc Oortl . Il968f) . This rapid 
rearrangement of photons in frequency space does 
not occur for photons that are far in the wing of 
the line profile, for which |a;| 3> 1. Here, the prob- 
ability that a photon is scattered from frequency x 
back into the line co re in a sin gle scattering event is 
extremely small fe.g. lHummeri . il962il . Instead, scat- 
tering pushes the photons back to th e line center 
by an average amount of only —1/x (jOsterbrock . 



(5) 



19621 ). Therefore, once a photon finds itself far in 
the line wing, it can stay there for ~ scattering 
events ( Adamd . 1972 ). It follows that for a popula- 
tion of Lya photons that are all in the blue wing of 
the line that: (i) < 0, and (ii) scattering does not 
immediately redistribute the photons in frequency 
space such that — s- Tgas. Therefore, a severely 
blueshifted Lya line has the potential to sustain a 
negative color temperature while it interacts with the 
gas. This negative color temperature may lead to a 
negative spin temperature, as we discuss next. 
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3.2. The Hydrogen Spin Temperature 



The expression for the equihbrium Hydrogen spin 



temperature is (see Appendix lA.ip 



^cmb Va'-^a H~ Vc^g. 



Vc 



(7) 



where = Va^ItI and = a^^,- Here, 
A21 and B21 are the Einstein rate coefficients and 
C21 is the colhsional de-excitation rate. This ex- 
pression differs shghtly from the standard expres- 



sion fe.g. iMadau et al.l . 119971 ). Here, Q[Ja(y)] is a 



function that depends on the exact Lya spectral 
shape (see Eg. IA.15"|) . and which classicaUy obtains 
a value Q[Ja{v)\ ~ 1. However, Hirata (2006, his 
Eq. B.18) has recently shown that due to a quantum 
interference term, the probability that a spin-flip 
occurs through scattering by Lya wing photons is 
strongly reduced. Therefore, Q[Ja{v)\ may become 
very small in cases in which all Lya photons are in 
the wings of the line profile. In this paper we will 
express our answers as a function of Q[Jq,(z^)], and 
keep in mind that it may obtain very small values. 

The denominator of Eq. [7] shows that > 
27/4(r„/r»)A2i - 3 X 10-11 s-i(T„/hioo 
K]){1 /Q[Ja{i^)]) results in T^ < 0, provided that 
the Lya scattering rate dominates the collision rate 
and the CMB scattering rate. The Lya scattering 
rate may be written as 



J{x)(j){x)dx, 



(8) 



where J^{x) is the normalised spectrum (J J'{x)dx = 
1). We will use this equation to get a constraint on 
the Lya luminosity in § 



4. The Lya Color Temperature in a Neutral 
Collapsing Gas Cloud 

We use th e Mo nte-Carlo code described in 
Diikstra et al. for our radiative transfer cal- 

culation. This code calculates the Lya transfer 
through spherically symmetric gas clouds, with ar- 
bitrary density and velocity fields (for a detailed 
description the interested reader is referred to Di- 
jkstra et al 2006). For the work presented in this 
paper the code was expanded to follow the time 
evolution of the Lya radiation field as a function 
of position. This was achieved by keeping track 
of the total time that elapsed since a Lya photon 
was emitted. When this time reaches user-specified 



values, a 'snapshot' is taken and the locations and 
frequency distribution of Lya photons in the frame 
of the inf ailing gas are recorded. 

4.1. The Model 

We study a Lya source point source that is em- 
bedded in a neutral collapsing gas cloud. The total 
(dark matter + baryons) mass of the halo hosting 
the collapsing cloud is Aftot = 10^^ Mq. Collapse 
occurs at redshift z = 8. Within the virial radius, 
Tvir = 16 kpc the gas density is assumed to trace 
the dark matter (which is given by a NEW profile 
with concentration param eter c = 5) modified by 
a core at r < 3rs/4 (see iMaUer fc BuUockL 12004 . 
their Eq. 9-11) where rg — ^vir/c. The infall veloc- 
ity field, w(r), is assumed to be linear with v(r) = 
^circ [r/ryir] in which Wcirc(=158 km s ^) is the cir- 
cular velocity of the host dark matter halo. This ve- 
locity field corresponds to that of a spherical top- 
hat model ( Diikstra et al.l . [2006h . For realistic initial 
density profiles, the mean density enclosed within 
radius r decreases smoothly with radius, causing the 
inner shells to be decelerated more relative to the 
overall expansion of the universe. This enhances the 
infall speed at small r relative to that in the top-hat 
model, and results in a fiatter velocity profile, pos- 
sible with an increasing infall speed towards smaller 
radii. Such a velocity field would result in a larger 
blue-shift of the Lya line (see Fig 4 of Dijkstra et al 
2006), which would not affect our result regarding 
the sign of the Lya color temperature. Beyond the 
virial radius, r > rvir, infall of gas continues and we 
adopt the density and velocity profiles based on the 
curves shown in Figures 2 and 4 of Barkana (2004; 
see Dijkstra et al, 2007, for a more detailed discus- 
sion of this choice). 

The gas collapsing outside the galaxy is assumed 
to be cold, Tgas — 300 K. This low gas tempera- 
ture reflects the fact that gas in the pre-reionisation 
universe has a low temperature . Furt hermore, cal- 
culations by iBirnboim fc Dekell (2003) have shown 
that virial shocks are likely not to exist in objects 
with masses below Mtot '--^ 2 x lO^^AfQ (also see 



Keres et al.l . 120051 ) . In the absence of a virial shock 



gas accretion occurs 'cold'. Of course, density in- 
homogeneities in supersonically collapsing gas will 
likely cause it to shock heat to a higher tempera- 
tures. However, instead of discussing whether our 
model is realistic or not, we simply point out that 
the Lya cross section is independent of gas temper- 
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Fig. 1. Time evolution of the hya spectrum 'seen' by a 
gas element at r = 15 kpc, inside a collapsing cloud of 
neutral hydrogen. The Lyo radiative transfer was calculated 
using a Monte-Carlo code. All Lya photons are inserted at 
t = 0, and the spectra are shown at t = 4fcross (black) and 
t = lOtcross (red) , where tcross = rsn/c is the light-crossing 
time to radius ran = 15 kpc. All spectra are systematically 
blueshifted relative to the line center, which is due to energy 
transfer from the infalling gas to the hya photons. These 
distributions of Lya photons have ^ —160 K (t = 4tcross) 
and Ta ~ —230 K (t = lOtcross)- When the hya scattering 
rate exceeds, Pa ^ 5 — 7 X 10~^^ the spin temperature 
may be negative, Ts < 0. 

ature in the wings of the hne profile. Therefore, the 
scattering process and the Lya color temperature 
are practically independent of the assumed gas tem- 
perature throughout this paper. 

Photons are injected at t = following a Gaus- 
sian frequency dist ribution with (t„ — Wcirc (see 
Diikstra et all 120071 for a motivation of this choice). 
Note that our results do not depend on the precise 
input spectrum of the Lya photons, as long as the 
collapsing cloud is optically thick to all emitted 
photons. 



4.2. The Results 

Snapshots were taken of the Lya frequency dis- 
tribution at Tsn = 5, 10, 15 and 20 kpc at four 
different times: t = 4icross, t = 6icross, t = Sicross 
and t — lOicross, where Across = ^sn/c is the light- 
crossing time to radius Vgn- The spectra 'seen' by 
the gas at r = 15 kpc are shown in Figure [1] as the 
black (t = 4tcross) and red {t = lOtcross) histogram. 
The Lya photons lie systematically to the blue side 
of the line. The inset of the Figure zooms on the 
line core and shows more clearly that no photons 



are present in the core of the line profile. The reason 
for the large shift to the blue side of the line can be 
easily understood: as the Lya photons propagate 
outward through an optically thick collapsing gas 
cloud, energy is transferred from the gas to the pho- 
tons, which results in an overall blueshift of the line 



2002t iDiikstra et al 



( Zheng fc Miralda-Escude . 
20061 ). This blueshift can exceed the infall velocity 
of the infalling gas significantly (large velocity shifts 
are common in Lya radiative transfer through op- 
tically thick gas clouds. In order for Lya to escape 
from an optically very thick gas cloud, the photons 
need to diffuse far into the line wings, where the gas 
cloud becomes optically thinner). 

At each subsequent timestep the Lya line shifts 
slightly further to the blue and broadens. This 
broadening of the line is a direct consequence of the 
frequency diffusion associated with the Lya scat- 
tering process. Frequency diffusion in a collapsing 
gas cloud preferentially shifts Lya photons further 
to the blue. As we discuss next, the Lya color 
temperature is negative at all times. 



4.3. The Lya Color Temperature and Hydrogen 
Spin Temperature 

Using Eq © we find that T„ = -160 K at t = 
4<cross and Tq = —230 K at t = lOtcross for a gas 
shell at r = 15 kpc (Fig. [1]). Therefore, neutral col- 
lapsing gas clouds surrounding a central Lya source 
may 'see' a negative Lya color temperature. At later 
times {t > lOicross) the Lya scattering rate de- 
creases to zero, as all photons have diffused away 
spatially. For a source that emits Lya photons at a 
constant rate, a gas element at r = 15 kpc would 
continuously be exposed by a weighted average of 
the spectra shown in Figure [TJ which would clearly 
also have an associated negative color temperature. 

For Ta = —200 K the spin temperature is neg- 
ative when the Lya scattering rate exceeds Pa ^ 
5 - 7 X 10-^^ /Q[Ja{iy)] s^i (§[321). This requires 
the Lya luminosity of the central source to exceed 
a few times W^^ / Q[Jaiiy)] erg s'^ (Eq.©. To pro- 
duce Lya luminosities of this magnitude requires 
high star formation rates. If stars form according 
to the conventional Salpeter mass function, the 
Lya luminosity and star formation rate, M*, are 
related by Llvq 1 - 3 xlO^ {M^/Mq yr ^) erg 



(e.g. I Diikstra et all . l2007l ). For a top-heavy 



mass function of stars, the Lya lu minosity may b e 
higher by an order of magnitude (|Schaereil[2003h . 
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Therefore, LLya ^ 10"^^ / QlJaii')] erg requires 
> 10'^/Q[Ja{iy)] Mq/jt. Note the importance 
of quantum interference for this work. As is shown 
in Appendix lA. 3) existing calculations suggest that 
Q[Jq(i/)] ~ 10^^, which implies that unphysically 
large star formation rates that are a factor of ~ 10^ 
times larger are required for a population inversion 
when quantum interference is included. 

A caveat to the previous discussion is that 
for star formation rates mentioned above, Af, > 
10'^/(5[Ja(z^)] M(7)/yr, the surrounding gas is very 
unlikely neutral and cold (which was assumed in the 
radiative transfer calculation) . This implies that al- 
though the scenario described above may result in a 
negative Lya color temperature, the Lya luminosity 
that is required to couple the Lya color tempera- 
ture to the hydrogen spin temperature is too large. 
In the next sections, we discuss additional factors 
that determine whether 21-cm masers can exist. 

5. The Radiation Pressure Exerted by 
Luminous Lya Sources 

We showed above that Lya luminosities exceed- 
ing ~ 10^^/(5[Jq(i^)] erg are required for a pop- 
ulation inversion in the hyperfine levels of atomic 
hydrogen in the collapsing gas. Here we show that 
the radiation pressure exerted by the Lya photons is 
likely large enough to halt the collapse of the cloud. 
The condition for Lya radiation pressure to unbind 
gravi tationally collapsed gas is fe.g .Oh fc Haiman. . 
l2002h 

GM^{<r) 

^LvQ^trap > C^bind ~ , (9) 

r 

where C/bind is the gas' gravitational binding en- 
ergy, M{< r) is the total mass enclosed in a 
sphere of radius r, and ttrap is the total time 
Lya photons are 'trapped' inside the sphere of 
radius r, before escaping. If the gas were opti- 
cally thin to Lya, then itrap = In optically 
thick gas, scattering of Lya photons increases this 
timescale to ttrap = /(''') x (^/c), where /(r) ^ 
15(T(r)/106)i/3(Tgas/10'* K)i/6, wh ere r(r) is the 
total Lya optical depth to radius r (|0h &: Haimanl . 
|2002[ ). For the density profile used in the radiative 
transfer calculation (§ 14. 3p we find that radiation 
pressure can exceed gravity for Lya luminosities 
Lhya ^ 10^^ erg s~^ (at all radii of interest). It is 
likely that while the central source's Lya luminosity 
is building up to the values required to establish a 
population inversion, the collapse of the cloud has 



halted. Since gas collapse was essential in the mech- 
anism proposed here to make Tq < 0, we conclude 
that the required Lya luminosities are too high 
(especially for very small values of QlJaii^)], which 
are favored by the most up-to-date calculations, see 



6. Constraints on the Properties of the 
Central Lyct Source 

A key requirement for a negative hydrogen spin 
temperature is the absence of -or more accurately, a 
low flux of- Lya photons in the line center, as these 
would scatter and rapidly rearrange the level pop- 
ulations in the hyperfine transition such that -—^ 
Tgas. So far, our calculation assumed that all Lya 
photons were emitted by the central source. Energy 
transfer from the gas to the photons caused Lya 
photons to blueshift into the line wing on their way 
out of the collapsing gas cloud, causing a large vol- 
ume of the cloud to practically see only Lya wing 
photons. However, in astrophysically realistic situa- 
tions Lya core photons may be injected at any loca- 
tion through different channels. These channels are 
discussed below. 



6.1. Alternative Sources of Lya Core Photons 

Different channels that produce Lya core photons 
include: 

1. Lya that is produced following photo- 
excitation by soft UV-photons, defined as photons 
in the energy range E ^ 10.2 - 13.6 eV. Soft-UV 
photons can generally penetrate deep into a neu- 
tral medium, where they are absorbed by neutral 
hydrogen atoms. Recombination events produce 
higher-order Lyman-series photons, which have a 
~ 35%/5 = 7% probability of being converted into 
a Lya core photon, per scattering event. The factor 
of 5 reflects the fact that a higher order Lyman- 
series photon scatters on average ~ 5 times before 
undergoing a cascade (which results in a Lya pho- 
ton ~ 35% of the time, see Pritchard & Furlanetto 
2006). 

2. Lya that is produced following x-ray heating. 
When hydrogen recombines after being photoionised 
by x-rays that penetrated the neutral gas, a Lya 
photon is produced with a probability of ~ 2/3 (for 
case B recombination). The emergent Lya photon 
has a high probability of being injected close to the 
line center. This effect is probably not important 
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because the recombination time tree — '^jineCt) ^ 
Across, particularly in an almost completely neutral 
medium (regardless of the gas temperature) . A more 
important effect however, is that photoionisation by 
x-rays may produce energetic electrons that excite 
the n = 2 level of atomic hydrogen, which immedi- 
ately resultsjii_the_jni^ core photon 
fe.g. lChen fc Miralda-EscudS l2006l) . 

3. Lya produced in collisions that excite the Lya 
transition. This effect is only important when the 
gas temperature is close to 10"' K. 

The newly produced Lya-core photons inter- 
act strongly with the neutral gas, and will drive 
th e spin t empera ture to th e gas temperature (as 
in IPeguchi fc Watsonl . 1985h . The evolution of the 



hydrogen spin temperature ultimately depends on 
which population of Lya photons dominates the 
spin-flip rate. If the central source is sufficiently 
bright in Lya relative to the soft-UV and x-ray flux, 
then the rate at which Lya wing photons produce 
spin-flips can still dominate. The prominence of the 
Lya emitted by a source relative to its continuum 
is quantified by its equivalent width. Therefore, the 
requirements mentioned above can be quantified as 
a constraint on the source's Lya equivalent width. 

6.2. Constraints on the Lya Source's Lya 
Equivalent Width 

The scattering rate of the newly created Lya core 
photons at radius r can be estimated as 



Pa.core(r) - {N)e{r), 



(10) 



Here, e(r) is the rate at which Lya core photons are 
injected at radius r. Furthermore, (N) is the total 
number of times a Lya core photon scatters before 
it disappears. To estimate (N), we first note that it 
takes on average ~ 10"' — 10^ scattering events for 
a Lya core ph oton to b e scatt ered into the wing, 
|a;| > xt 



3 ( Hummei 



19621 ). However, photons 
close to the core- wing transition frequency, \xt\, can 
scatter back into the line core in a single scatter- 
ing event. Therefore, (N) should be higher. Veloc- 
ity gradients facilitate the escape of Lya photons. 
In the presence of velocity gradients the majority of 
the scattering events of a photon occur around its 
resonant region. For a static medium, a Lya pho- 
ton scatters on av erage ^ tq times as it traverses a n 
optical depth tq ( Adamd . 19721 : Harrington . 1973 ). 
The total optical depth of the resonant region, Tres, 
can be calculated from the total optical depth over 



the scale for which the velocity changes by a ther- 
mal width (also known as the Sobolev length), vth- 
Tics = nHcroVth[dvin{i,ii/dr]~''- . If we approximate 



z/fvir, where Ucirc and Tvir are the 



rfWinfall/rfr ~ 

circular velocity and virial radius of the host dark 
matter halo, then r^os ~ 2 x 10*^ [(1 + (S)/60][(l + 
Zvir)/9]'^/^. Here, 5 is the overdensity of the gas. 
Therefore, for a photon to escape from its local res- 
onant region, it scatters on average {N) ~ Tj-cs ~ 
10^ - 10'^ times. 

In § 14.31 we showed that < required Pa ^ 
5 — 7 X /QlJaiv)] s~^ Let us require that the 

scattering rate of Lya core photons is subdominant, 
e-g- -Pq, coro('") < 10"^^ s~^ (the factor Q[Ja{v)\ is 
deliberately omitted, as quantum interference in the 
core can be ignored) . Then, Eq llOl implies that e(r) < 
10-" s-V(7V). 

The rate at which soft-UV photons are converted 
into Lya is given by 



e(r) - 0.07 x Att 



Ly7 



hv 



■a,,t{v)e-^^^^^~>dv. (11) 



Here, the factor 0.07 is the mean probability for a 
higher-order Lym an-series photon to be c o nverte d 
per scattering ( Pritchard fc Furlanettol . l2006l ). 
t(i/, r) = NH{r)(Ttot{v), in which crtot(j^) is the 
absorption cross section at frequency v which is ob- 
tained by summing over all Lyman series transitions 
(derived in Appendix [B| , and Nh {r) is the total 
column of neutral hydrogen between the source 
and radius r. The integral is performed over hv G 
[12.8 — 13.6] eV, as Ly/3 photons cannot be con- 
verte d into Lya photons ( Pritchard fc Furlanettol 
l2006f ). We assume that the flux density is con- 
stant, i.e. J{v) ~ K, where i^T is a constant, and 
that Nnir) — 10^^ cm"^ and numerically eval- 
uate the integral. The constraint e(r) < 10^^^ 
s-V(A^) then translates to K < 1.5 x 10~'^'^/{N) 
ergs s"^ cm~^ sr^^ Hz^^, which can be recast as 
IGn'^r'^K < 10^^ /{N) ergs s^^ A^^. Since the Lya 
luminosity of the central source was required to be 
Lhya ^ 10"'^/Q[Ja(i')] ergs s~^, this requires that 
the Lya equivalent width is EW> O.l(iV) /Q[Jaiiy)] 
A= 106/Q[J„(j/)] A, for (N) = 10^. Although the 
precise lower limit depends on model parameters 
such as A^Hi and (A^), our lower limit on EW pro- 
vides a reasonable order-of-magnitude estimate. 

For comparison, a normal population can theoret- 
ically produce a maximum of EW=200-300A (e.g. 



Schaereij . 120031 ) . This number can be higher by about 



a factor of 5 for population HI stars and quasars, but 
this is still well below what is required to produce 
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a 21-cm maser. However, it has been shown that a 
muhiphase intersteUar medium (ISM) can preferen- 
tiaUy transmi t Lyg photons over contin uum pho- 
tons (|Neufeldl . [l99ltlHansen fc Ohl . l2006l ). In multi- 
phase ISM models, the ISM consists of cold, dense, 
dusty clumps embedded in a hot medium. The Lya 
scatters off the surface of the cold clumps and prop- 
agates mainly through the hot, dust free, interclump 
medium. On the other hand, the continuum photons 
penetrate the clumps where it is subject to absorp- 
tion by dust. Therefore, a multiphase ISM can pref- 
erentially absorb the continuum flux and boost the 
Lya EW. It is currently not known whether this ef- 
fect actually occurs in nature, especially at the level 
that is required to make EW> lQ^/Q[Ja{v)] A. We 
point out that the multi-phase ISM that is required 
to produce EW> IQ^ / Q[J a{v)] A is likely also suf- 
ficient to counteract the negative impact of x-rays. 
For a spectrum of the form cx u^"' the injection 
rate of Lya photons through x-rays is comparable 
to that given by Eq (fTTj) for a > 2. 



7. Conclusions 

We have performed Monte-Carlo calculations of 
the radiative transfer of Lya photons emitted by a 
source embedded in a neutral collapsing gas cloud 
and computed the Lya spectrum as function of 
radius and time. This source represents a young 
galaxy or quasar during the early stages of the epoch 
of reionisation (EoR). After deriving an expression 
for the Lya color temperature Tq for arbitrary fre- 
quency distributions of Lya photons (Eq. [H]), we 
show that neutral collapsing gas clouds surrounding 
a central Lya source may see a negative Lya color 
temperature (Fig[T]) under certain circumstances. 

We have investigated whether a population inver- 
sion in the hyperfine transition of atomic hydrogen 
via the Wouthuysen-Field effect can be established 
in such scenarios, thus producing a 21-cm maser. 
The reason for this investigation is clear: if 21-cm 
masers exist during the epoch of reionisation (EoR) , 
this could greatly boost the expected 21-cm emis- 
sivity of neutral gas during the EoR. Unfortunately, 
we have found that in practice a population inver- 
sion does not occur for several reasons: 

(1) The required Lya luminosities exceed ~ 10**^ 
erg s~^. Although this can be obtained with phys- 
ically possible star formation rates (M, > 1{)^Mq 
yr~^) the radiation pressure exerted by such a large 
Lya flux likely halts the collapse of the cloud (§ [5|) . 



(2) Quantum interference can strongly reduce the 
strength of th e WF-effect in the Lya line wings 
()Hiratal . 2006f ). The probability that a spin- flip oc- 
curs when a Lya wing photon scatters is reduced 
by a factor of Q = {x^^/xY (see Appendix IA.3P 
when quantum interference is included. In the con- 
text of this paper, this factor can be as small as 
Q = 10^"* — 10~^. Therefore, when this quantum 
correction to the WF-coupling strength is applied, 
the required Lya luminosities are larger by orders 
of magnitude. 

(3) A key requirement for having < is a very 
low level of Lya photons in the line center, as these 
would scatter and rapidly rearrange the level pop- 
ulations in the hyperfine transition such that 
Tgas- We discussed in § [6l however, that Lya core 
photons can be injected at a given location via differ- 
ent channels including x-ray heating, coUisional ex- 
citation and conversion of higher Lyman series pho- 
tons. The negative impact of these newly produced 
Lya photons (both by x-rays and soft-UV photons) 
is only overcome when the central source has an ex- 
tremely prominent Lya line with EW> lO^/Q A, 
which is several orders of magnitude larger than that 
of currently known sources in the Universe. 

Unless the high redshift Universe contains some 
truly remarkable sources that can overcome the 
three obstacles mentioned above simultaneously, 
the Wouthuysen-Field effect will not produce a pop- 
ulation inversion in the 21-cm transition. Detecting 
the redshifted 21-cm from the epoch of reionisation 
will therefore remain among the greatest observa- 
tional challenges in cosmology. 
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Fig. A.l. Schematic depiction of hyperfine splitting of the electronic ground state (n=l) and the first excited state (n=2) of 
the hydrogen atom. This 6-leveI model has been used to arrive at the expressions for the Lya color temperature, Eg. lA.Tl and 
the spin temperature, Eg. IA.9I The dashed lines denote the transitions involved in the mixing of hyperfine levels. Levels 3 and 
6 do not contribute to the mixing of levels. The notation used for the individual levels is pLj, where L is the orbital angular 
momentum, J the total angular momentum quantum number. Here, F is the sum of J and I, the nuclear spin quantum number. 



Appendix A. Derivation of the Lya Color and Hydrogen Spin Temperature 



A.l. The 6-Level Hydrogen Atom without Collisions and CMB photons 



In this section we derive the expression for the Lya color temperature taking into account hyperfine splitting 
in the Hydrogen atom's n = 1 and n = 2 states, following Meiksin (2000). Ignoring collisions and absorption 
and re-emission of CMB photons, the rate equation for atoms in the n = 1 state is 



dfii 



J 



B 



A 



42 



14- 



fli + B2 



A, 



= ^2^21 + J 



A42 + A41 

54 A41A42 



712 - Bi 



A 



52 



-ni 



gi Ai2 + Ail 52 -441 + A. 



Ail + ^42 

^ gi Ai2Aii 

42 



-n2 



A52 + ^51 
55 A51A52 



rii + B; 



A^i 



25 



n2 



-ni 



A51 + ^52 
55 A52A51 



gi A52 + A51 52 ^51 + ^52 



"2 



(A.l) 



where we used the relations B^iu = j^B^d = fj^fi^^"'^' ^'^'^ notation J = J J{v)(j){v)dv. Following 



MeiksinI (|2000f) . 



Ail — g^tt' ^42 



-:A„ 



A51 = A52 



(A.2) 



where A^ — 6.25 x 10^ s is the total Einstein A-coefficient, and gi = 1, 52 = 3, 54 = 3 and g^ = 3. We get 



drii 



= »^2^21 



2hv^ 



Aa X 



(A.3) 
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-ni / J{i')(l)ii{i/)di' + -n2 / J{i^)(f>2iii^)dh' — -ni / J{v)4)iz{v)dv + -n2 / J{v)4>25{v)dv 
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All line profiles are identical, apart from an offset in frequency: (j)2i{i^) — <^i4(' ^ + At/ . ^), (jj^h li^) — </'i5(i' + A^^g) 
and (l)ii{iy) = (t>i5{v + Ave), where Ai/q — 10.89 Ghz and Az/5 = 1.42 Ghz ( Hiratal . [2006l ). If we approximate 
(f>{h' + Aiy) w 0(i^) + Av^ and write everything in terms of (which is written as ipi^) hereafter), and 
replace n2 with uh — ni, we get 



drii 
'dT 



n2A 



21 



9hi 



>Aa X 



(A.4) 



UH J JH(2<Piiy) + ^AiyG + 2^A,,s)dt^~ni J J{v) 
This can be written in the following more compact form 
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{v) + ^^AvG + 2-^ Ai^s dv 



^ = (1 - y)A2i + 61 - 62 - y(46i - 62) (A.5) 

after dividing the left and right hand side of Eq lA.Sl bv nn and after writing y ~ ni/nn and where we defined 

f / d(h \ 2c^ f 86 

6r - j J(.)(20(.) + ^A.,)d.; ^ ^A.^A.s J £jMd. (A.6) 

The equihbrium value of y is obtained by setting dy/dt—0 and corresponds to a spin temperature Tg of 



3w 61 

4y - 1 62 



^ (l){v)J{v)dv AvG 



/ (t){v)J{v)dv 



(A.7) 



where we used thatT, < TsII]. Furthermore, we have used that the second term within the square brackets is 
negligible compared to the first term for all applications in this paper. Here the last definition only states that 
if Lyck scattering were the dominant process regulating the level populations in the 21 cm transition, then the 
equilibrium spin temperature is equal to the Lya color temperature. Since hAvs = ksT.,, this yields Eq. ([5]). 



A. 2. The 6-Level Hydrogen Atom with Collisions and CMB photons 

In the presence of CMB photons and when collision are accounted for, Eq. (|A.5[) changes to 

^ = (1 - 2/) (^21 + C21 + S2i/cmb) + ibi - 62) - y (4^1 - 62 + C12 + Bi2/cmb) • (A.8) 
In this case the equilibrium solution (obtained by setting dy/dt=0) corresponds to a spin temperature T, of 



42/ -r 



A21 



C21T, 



biT, 



C21T, 



, (A.9) 



A2lTa A2lTgaB 



where we used that B12 — 3i?2i and C12 = 3C2ie ^*/'^8''= « 3C2i(l — T^/Tgas). When the Lya scattering 
rate is expressed as 



An 



hv 



(j{v)dv ; 



47r BahvQ 
hvQ 47r 



J [v)(j){v)dv ~ SAa 



27 

2K' / ^^-^"^^ 



the standard equation for the hydrogen spin temperature fe.g. iMadau et al.l . ll997f) is recovered: 



^* ~l~ ^cmb ~t~ UoTcx ~l~ Uc'-^gi 

i + ya + yc 



27^2 iTa 



C'2iT', 

^2l7Ka; 



(A.ll) 



where in the numerator is usually omitted as it is negligibly small compared to the other terms in most 
applications. 



A. 3. Expressions ofT^ andTg with Quantum Interference 



Hirata ( 20061 ) has performed detailed quantum mechanical calculations of the cross-section of the process 



in which absorption and re-emission of Lya photons results in a transition between levels 1 and 2 (Fig I A. l]) 



^ At first glance, it is not obvious that this inequality is satisfied. The spin temperature is usually Ts 2> T«, while in this 
paper Ts <S T*. One may therefore expect that there is a regime in which Ts ~ 0. However, this is not the case. The main 
reason is that continuously increasing the number of atoms in the n = 2 state, say from a state in which Ts 

— ^ gas 1 increases 

Ts, where in the limit 712 — » 3nH/4 one gets Ts — > 00. Increasing the fraction of atoms in the n = 2 state further results in a 
change in the sign of the spin temperature and Ts — > —00. By further increasing the number of atoms in the n = 2 state, one 
reduces the absolute value of Ts. In practice, the limit Ts ~ is never reached as it would correspond to a case in which all 
atoms are in the n = 2 state, i.e. n2 — > nj]- 
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These cross-sections are denoted as ai2 and (T21 and represent the transitions 1 — > 2 and 2 — > 1, respectively. 
They are given by 

2 2 

Cri2 = gCTo[014 + (1)15 - 2015/14], 0-21 = 27^0 ["^24 + 025 ^ 2025/24]- (A. 12) 

Here, (Jq is the Lya cross-section evaluated at resonance and 0i5/i4 and 025/24 are quantum interference 
term££] which are defined as 

0i5/i4('^) ^ /TT- due ^ ''0i5/i4(i^ ); 0i5/i4('^) = -T TT T' 



(and 025/24 is defined equivalently) where ctv — \J ^ y ^o, in which mp = 1.66 x 10 ^'^ g is the proton mass. 
Further i nore, 7 = Aa/'ii:, Ai/15 = (y — 1^15), where vi^ is the frequency at resonance of the 1^5 transition 
(jHiratal . [2OO6I . ;he other symbols are defined simi larly). For comparison, a no rmal 'unconvolved line profile' 



is given by a Lorentzian 0i5(i^) — a i , 2 (e.g. iRvbicki k, LightmanL 19791 . Eq 10.73). For Lya scattering 



in the line core, the quantum interference term is typically smaller than the traditional terms and it can be 
ignored without introducing a large error. However, this is not true in the far line wings. Here, by combining 
Eq |A.12l and Eq lA. 131 the line profiles 0i2(a;) and 02i(a;) can be written as 

, 4a/ 0:45 \ 2 A a ( X45 \ 2 ■ r / a , < \ 

012 = - ^ 2 1 j ' 012 = ^ 2 1 j ' with quantum interference, (A. 14) 

where 0:45 = (z^i4 — v\^l Avu- For comparison 

4 a ,4a . , • . 

without quantum interference. 



9^AFa;2' 27- 



Therefore, the probability that a spin-flip occurs when a Lya wing photon scatters is reduced by a factor of 
(x45/a:)2 when quantum interference is included. In this paper, this factor can be as small as 10"'' - IQ-^! We 
caution that the exact degree at which quantum interference reduces the spin-flip probability is uncertain. 
The main reason for this uncertainty is that the unconvolved line profiles may deviate from Lorentzian in the 
far line wings, which may affect these calculations (Zygelman, private communication). 

Regardless of whether quantum interference is included, we can write 021 (i^) = <yi\{v)loQ — i0i2(i^ + Ai/g). 
By repeating the analysis of i jA.ll one can show that the expression for the Lya color remains the same. 
However, in the expression for the hydrogen spin temperature ya yaQ[Jct(j^)] where 

ntT / \^ / ■^a(t^)[024 + 025 - 2 025/24] rfi^ 
J JaWjmi + (P25\dv 

This equation simply states that the cross-section that includes quantum interference should replace the 
cross-section that ignores this term in the calculation of the spin-flip rate. Note that if all Lya were far in the 
line wings, then Q[Ja] oc x~'^ . 

Appendix B. The Opacity of Gas to Photons in the Energy Range E — 10.2 — 13.6 eV. 

The cross section for Lya absorption can be written as 

7re2 1 

cTa{x) = cra,o0(a;;yl2i) = /12 — 0(a;; A21), (B.l) 

mc -y/TrAi/o 

where an n is the line center cro ss section, /12 = 0.4167 the Lya oscillator strength and (j){x; A21) is the Voigt 



function ( Rybicki fc LightmanL [l979i . p. 288, Eq. 10.70). The reason for writing the voigt function as 0(a;; A21) 



^ Note that when the quantum interference terms are ignored the rate equation for ni can be written as = n2{A2i + 
P21) — niPi2, where P12/21 is the rate at which Ly« scattering puts atoms from the 1 ^ 2/2 — > 1 state. This can be written 
as^^ = n2(^2l +iTvn2 J '^j^"J 021 {v)dv — 47rni J' '^j:^J 021 {u)du. If we use Ea lA.12l and ignore the quantum interference term, 
we are left with E g IA.4I Therefo re, apart from the quantum interference term, our calculations are in exact agreement with 
those presented by iHiratal 1 I2OO6I) . 
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Fig. B.l. The absorption cross section of neutral atomic hydrogen gas to soft UV photons. The spikes at E=12.1, 12.8, 
13.05, ... correspond to the Ly/3, Ly7, LyiS,... -resonances, respectively. The horizontal dashed line corresponds to the Thomson 
cross-section. We assumed that Tgas = 300 K, but note that the gas temperature barely affects this plot (see text). 

is that it explicitly demonstrates that it depends on the Einstein A-coefBcient (in the wing of the line profile, 
(p{x) = al\\pKx^\, where a — Ai^j ^-k^vj)). Similarly, the cross section for higher order Lyman-series photons 
can be expressed as: 



^n^'^ri) — fln~ 



1 



(B.2) 



mc y^AVn 

Here, Xn = {v — Vn)//^.^,^. To compute the cro ss sections therefo r e req uires knowledge of the oscillator 
strengths, /i n, and the Einstein coefficients, A„i. Menzel fc Pekerij ( 1935 ) derived the oscillators strengths 
(Eq. 10.46 in lRvbicki fc Lightman, .1979.1 to be 

_ 29n5(n- l)^"-4 

" 3(n-M)2"+4 • ^^-^^ 
From the definition of the oscillator strength we know that 

Tre^ hVn gn hVn „ 9n hVn . (J. 
Jin — —. i'lri — -. i>jil — ^nlTTJ — T 

mc An gi An gi An ^hVn 

(jRvbicki fc Lightmanl . ll979l . Eq. 10.70). A combination of Eqs. (|B.3p and (|B.4p yields v4i„. 

Given the cross-sections for the various Lyman-series transitions, it is possible to calculate the opacity of 
neutral gas for a photon of any given energy in the range E = 10.2 — 13.6 eV, by summing over all Lyman- 
series transitions 



i=l 



a.i{xi), Xi = {v - Vi)//^Vi 



(B.5) 



The result of this calculation is shown in Figure [BTT] Here, we assumed that Tgas = 300 K. However, the 
appearance of this plot is very insensitive to the gas temperature, as this only affects the heights and widths 
of the resonance spikes in very narrow frequency intervals. 
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